Reduction of the hydrophobic attraction between charged solutes in water 



m 
o 
o 

(N 

» 

o 
O 
o 



J. DzubiellfEl and J. -P. Hansen 
University Chemical Laboratory, Lensfield Road, Cambridge CB2 1EW, United Kingdom 

(Dated: February 2, 2008) 

We examine the effective force between two nanometer scale solutes in water by Molecular Dynam- 
ics simulations. Macroscopic considerations predict a strong reduction of the hydrophobic attraction 
between solutes when the latter are charged. This is confirmed by the simulations which point to 
a surprising constancy of the effective force between oppositely charged solutes at contact, while 
like charged solutes lead to significantly different behavior between positive and negative pairs. 
The latter exhibit the phenomenon of "like-charge attraction" previously observed in some colloidal 
dispersions. 

PACS numbers: 82.70.Uv,87.16.Ac,61.20Ja,68.08.Bc 
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When apolar solutes are dispersed in water, they exert 
effective, solvent induced forces on each other, generally 
referred to as hydrophobic attraction. For large solutes, 
say globular proteins, the mechanism for hydrophobic in- 
teractions can be traced back to solvent depletion (or 
"drying" ) from the volume bounded by opposite surfaces 
of the solutes Q, 0. This mechanism, which is remi- 
niscent of polymer-induced depletion attraction between 
colloidal particles on larger (mesoscopic) scales [3j, is 
most pronounced when the solvent is near liquid-vapour 
coexistence 0. However, many biomolecular solutes are 
charged, and in this work we present evidence, based on 
extensive Molecular Dynamics (MD) simulations, that 
the electric field due to the charges carried by the solutes 
leads to a considerable reduction of the hydrophobic at- 
traction. 

Such a reduction may be qualitatively understood from 
macroscopic considerations. It is a well-known fact that 
a polar liquid, like water, will rise inside a vertical con- 
denser to minimize the overall electrostatic energy (the 
rise being limited by gravity on a macroscopic scale) Q . 
Transposing this situation to the nanometric scale, con- 
sider two parallel plate-like solutes of area A, separated 
by a distance D, inside a polar solvent of relative di- 
electric permittivity e, carrying opposite surface charges 
±<r. Neglecting edge effects, the electric field between the 
plates is Eq/c with E — cr/2e . We require the differ- 
ence in the grand potential between the situations where 
the liquid solvent (/) or its vapour (g) fill the volume AD 
between the two plates: 



saturation, one easily arrives at the following expression 
for the difference in grand potentials per unit area: 
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where P a is the pressure of phase a and j wa the surface 
tension between phase a and the plate ( "wall" ). Consider 
a state close to phase coexistence at temperature T, and 
let Sp — p — |Lt sa t be the positive deviation of the chem- 
ical potential from its saturation value. Expanding the 
P a to linear order in 5p around their common value at 
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At the "drying" transition between the plates, fli — £l g — 
and 7„i - j wg = ji g = 7. Since p g < pi, and e; = e > 
e g ~ 1, eq. (2) yields the following expression for the 
critical distance D c between the plates at which drying 
occurs: 
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For strong electric fields, (Eq < 10 10 V/m corresponding 
to surface charges a < e/nm 2 ), the electrostatic term in 
the denominator is typically 10 times larger than the bp 
term in the vicinity of gas-liquid coexistence, and leads 
to values of D c of the order of a few A. This strong re- 
duction of D c hints at a considerable weakening of the 
hydrophobic interaction between two solutes when the 
latter are charged. Note that within our macroscopic 
model this reduction is due to the overall electric field, 
not to any Hydrogen-bonding of the solvent molecules to 
hydrophilic "patches" on the solute surface. 

In order to confirm the qualitative prediction of the 
schematic model, we have carried out extensive MD sim- 
ulations of two spherical solutes immersed in a bath of 
SPC/E water molecules p|. The solutes are spheres of 
radius R which repel the solvent molecules by a repulsive 
e(r — R)~ 12 , where r is the distance from the solute center 
to the oxygen atom of a water molecule; the energy scale 
e is chosen such that the O atom experiences an energy 
k^T at a distance of r — R =lA from the solute surface. 
The simulation cell is a cube of length up to L = 40A, 
containing up to 2000 water molecules, depending on the 
solute size; solutes are placed at fixed positions on the 
body diagonal of the simulation cell. The box dimen- 
sions are chosen such that the surface to surface distance 
to the nearest image solute is at least 20A. 

The MD simulations are carried out with the 
DLPOLY2 package @, using the Verlet algorithm 0, 
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FIG. 1: Density profiles of the water molecules around 
two neutral spherical solutes of radius R = 10A (a) 
and two oppositely charged spherical solutes of radius 
R = 10A carrying a charge ±qe with (b) q = 2, 
(c) q = 5, and (d) q = 10. The surface-to-surface 
distance in all cases is s — 4A. In the contour plots 
dark regions show low density regions while high den- 
sities are plotted bright. The panels below the con- 
tour plots show the water density p scaled with wa- 
ter bulk density po in a cylinder of radius R c = 5A, 
coaxial with the center-to-center line of the solutes. 



with a timestep of 2fs. The Berendsen barostat and ther- 
mostat || were used to maintain the SPC/E water at a 
pressure of 1 bar and a temperature T — 300K. All elec- 
trostatic interactions were calculated using particle-mesh 
Ewald summations 0. We first consider the case of un- 
charged solutes. The water density profile is illustrated 
in frame (a) of Fig. 1 for the case of solutes of radius 
R = 10A and a surface-to-surface distance along the 
z-axis joining the centers s = 4A. The upper part of 
the frame shows a density contour plot coded by vari- 
able shades of grey. The lower part shows density pro- 
files along the center-to-center to axis z, averaged over a 
coaxial-cylindrical volume of radius 5A. The density pro- 
files show a considerable depletion of the solvent within 
a radial distance of 5A from the center-to-center axis, 
reminiscent of the observations of Wallquist and Berne 
for flatter solutes As the surface-to-surface distance 
s is increased for fixed radius i?, the water molecules pen- 
etrate into the region between opposite solute surfaces, as 
signalled by a rapid increase of the central peak (around 




FIG. 2: Simulation results (symbols) of the mean force, 
F* = /SPA, between two neutral (q = 0) spheres 
in SPC/E water. Error bars are omitted for clar- 
ity; the error is estimated to be AF* = 0.4, com- 
parable to the symbol size. The lines are guide to 
the eye. Results are plotted for solute radii R — 
3 A (circles), R — 5 A (squares), R — 8A (dia- 
monds), R — 10 A (triangles pointing up), and R — 
12A (t riangles pointing left). The inset shows the inte- 
grated force (potential of mean force) in obvious order. 



z = 0) in the density profiles. When s w 6.5A, the sol- 
vent layers around an isolated solute are hardly disturbed 
by the presence of the other solute. 

The mean effective force acting on each of the solutes 
in the presence of the second at a surface-to-surface dis- 
tance s is calculated by averaging the total force due to 
all solvent molecules over the configurations generated by 
MD runs extending over typically Ins. This average force 
obviously goes to zero at large distances s and for sym- 
metry reasons, it is directed along the center-to-center 
axis. Examples for several radii 3A < R < 12A are 
shown in Fig. 2. The largest radii are of the order of the 
size of small globular proteins or of oil-in- water micelles. 
As expected from a depletion mechanism, the force is at- 
tractive and its contact values and range increase with 
R. The potentials of mean force w(s) may be calculated 
for each R by integrating the force. The resulting po- 
tentials are shown in the inset to Fig. 2. They closely 
resemble results obtained for polymer-induced depletion 
potentials between spherical colloids, albeit on different 
length and energy scales 0, ^| . Note that the force at 
contact, F(0), scales roughly with R. This may be ra- 
tionalized by a simple consideration of the potential of 
mean force for plates (R — 00) at contact, w(0) = —27 
(where 7 is now the plate-solvent surface tension which 
differs little from the liquid-gas surface tension ^lj) and 
an application of the Derjaguin approximation, valid for 
weakly curved substrates (i.e. large R) 0]; this leads 
to the estimate F(s = 0) = — 2ttR'j, which indeed pre- 
dicts linear scaling; the value 7 ~ 0.05 J/m 2 , extracted 
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FIG. 3: Density profiles of water oxygen and hydro- 
gen atoms (inset) around one isolated solute with ra- 
dius R = 10 A and central charge q = (cir- 
cles), q = —10 (squares), and q = 10 (diamonds). 
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FIG. 4: Same as Fig. 2 but now the mean force, F* — 
/3-FA, is for oppositely charged spheres of R = 10A and 
different central charges ±qe are shown: q — (circles), 
q = 2 (squares), g = 5 (diamonds), q — 10 (triangles 
up). The dashed line represents the electrostatic force 
between 2 periodically repeated solutes with opposite 
charges q — ± 10 in a continuous solvent with permit- 
tivity e = 80. The inset shows the resulting potentials of 
mean force; the contact values w(s = 0) increase with q. 



from the MD value of F(Q) for the largest R, is reason- 
ably close to the liquid-vapour surface tension of water 
under normal conditions (7 = 0.073 J/m 2 ). The agree- 
ment must be considered quite satisfactory in view of the 
roughness of the estimate. 

We now turn to the charged solutes. Frames (b)-(d) in 
Fig. 1 show water density profiles in the vicinity of two 
spheres carrying opposite electric charges ±ge at their 
center (opposite charges ensure overall charge neutrality 
without any need for counterions). As q increases from 



zero (frame (a)), water is seen to penetrate between the 
two solutes, the central peak around z = in the den- 
sity profiles increases rapidly and its amplitude reaches 
roughly the bulk density of water when q — 10. Note 
that this central peak is asymmetrically split, indicating 
the presence of two hydration layers which differ some- 
what depending on their association with the anionic or 
cationic solute. This difference is also evident in the con- 
tact values of the outside surfaces of the solutes, and is 
a consequence of the different arrangements of the water 
dipoles around the solutes induced by the local electric 
fields. The asymmetry of the profiles can be rational- 
ized by inspecting the density profiles of O and H atoms 
around isolated solutes, plotted in Fig. 3. The hydration 
shell is more sharply defined around the cationic than 
around the anionic solute. The water dipoles tend to 
point radially away from the cation, while the opposite 
configuration is more favourable around anions. Note 
also that driving water into a narrowly confined region 
under the action of a strong electric field (here between 
oppositely charged solutes) is an effect reminiscent of that 
observed in recent MD simulations of ion permeation of 
hydrophobic nanopores |12j . 

The resulting mean forces between solutes are plot- 
ted for q = 0,2,5 and 10, as functions of the surface-to- 
surface distance s in Fig. 4 together with corresponding 
potentials of mean force. The mean force includes the di- 
rect Coulomb interaction between the two solutes (with 
proper account for the periodic images), which is in fact 
an order of magnitude larger that the total mean force. 
At large distances hydrophobic interactions become neg- 
ligible and the force should tend to — q 2 e 2 /(Aireoer 2 ), 
where r = 2R + s and e is the dielectric permittivity 
of bulk SPC/E water; the corresponding curve is also 
shown in Fig. 4. 

The most striking result illustrated in Fig. 4 is the 
near independence of the force at contact, s = 0, with 
respect to solute charge. From the density profiles in 
Fig. 1 the hydrophobic attraction is expected to be re- 
duced but this reduction is almost exactly compensated 
by the Coulomb attraction between solutes in the pres- 
ence of the solvent. As q increases, the initial slope of 
the effective force increases. The potential of mean force 
(shown in the inset of Fig. 4) exhibits a contact value 
which increases with q, indicating that the reduction of 
hydrophobic attractive free energy and the reduction of 
dielectric screening due to restructuring of water by the 
strong electric field of the solutes clearly outweigh the 
increase in bare Coulomb attraction between the latter. 
Simulations calculating the forces at and near contact 
for q — 7 and q — 15, not shown in Fig. 4, confirm this 
trend. Note that the potential of mean force curve for 
q = 10 shows more long range attraction as compared 
to the smaller q data due to the increased electrostatic 
attractive interaction. The eye-catching kink in the force 
for q =10, at a distance s rj lA is reproducible, and is 
probably related to the pronounced shell structure of wa- 
ter molecules around highly charged solutes, illustrated 
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FIG. 5: Same as Fig. 2 but now the mean force, 
F* = I3FA, is for different signs of the charged 
spheres of radius R = 10A and fixed central 
charges qe = ±5e are shown: oppositely charged 
(+ -, circles), equally charged with q = — 5 (- 
, squares), equally charged with q — 5 (++, diamonds). 



in Fig. 3. While for neutral (and weakly charged) so- 
lutes, the O and H density profiles show little structure, 
they are sharply peaked at a distance s « lA of the O 
atoms from the solute surface. This would lead to a com- 
plete shared hydration layer, and consequently to a kink 
in the force versus distance curve, between 2 flat solutes 
separated by s = 2A . This critical separation is shifted 
to shorter distances due to the curvature of spherical so- 
lutes. 

In view of this delicate balance between various in- 
teractions, we have also examined the case of equally 
charged solutes. In this case monovalent counterions 
(Na + or CI - ) were included to ensure overall charge neu- 
trality. The situation is summarized in Fig. 5 for solutes 



of radius R = 10A and charge q = ±5. The water density 
profiles (not shown here) are symmetric with respect to 
z = for equally charged solutes, but differ substantially 
when going from a pair of anions to a pair of cations. In 
the latter case water is much more structured into well 
defined hydration shells, as shown in Fig. 3. This dif- 
ference is reflected in the effective forces and potentials 
shown in Fig. 5. While the force and potential of mean 
force are practically zero at all distances in the case of 
cationic solutes, anionic solutes attract each other sig- 
nificantly showing that residual hydrophobic attraction 
overcomes the repulsion between like charges, which may 
be more efficiently reduced by the local permittivity of 
water in the immediate vicinity of the anionic solutes. 

To summarize, we have shown that electric charges car- 
ried by nanometer-scale solutes have a profound influence 
on hydrophobic interactions. Our MD simulations con- 
firm that for neutral solutes hydrophobic attraction is 
intimately linked to solvent depletion, and the contact 
value of the hydrophobic force is directly related to the 
surface tension of the pure solvent. Solvent depletion is 
suppressed by the electric field due to any charge carried 
by the solute, but for oppositely charged solutes the re- 
sulting loss of hydrophobic attraction is compensated by 
their mutual Coulomb attraction, leading to a nearly con- 
stant contact force. The effective force between equally 
charged solutes depends on the sign of their charge. An- 
ionic solutes exhibit a striking "like-charge attraction" , 
a subject of considerable debate in the field of meso- 
scopic charge-stabilized colloids Coulombic and 

hydrophobic effects are intimately entangled, and their 
inter-relation is expected to play an important role in the 
analysis of protein aggregation and related biomolccular 
mechanisms. 
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